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Abstract 

Contrarily to a common belief, any beam of light possesses to a some extent the ability to 
“reconstruct itself” after hitting an obstacle. The celebrated Arago spot phenomenon is nothing 
but a manifestation of this property. In this work we analyze the self-healing mechanism from both 
a mathematical and a physical point of view, eventually hnding a new expression for the minimum 
reconstruction distance, which is valid for any kind of beam, including Gaussian ones. Finally, a 
witness function that quantify the self-reconstruction capability of a beam is proposed and tested. 
The results presented here help clarifying the physics underlying self-healing mechanism in optical 
beams. 
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I. INTRODUCTION 


In this paper we illustrate, with the help of several examples, the self-healing property 
exhibited by different types of optical beams. In particular, we analyze and compare the 
behavior of Gaussian, Bessel, Bessel-Gauss beams and of a pair of plane waves. In recent 
years, the remarkable capacity of a beam to reconstruct itself after encountering an obstacle 
(self-healing mechanism) attracted considerable interest and has been the subject of numer¬ 
ous theoretical and experimental investigations HHZI. In this context, particularly relevant 
is the theoretical work by Ghu and coworkers Pi- 

The aim of this work is to study further the self-reconstruction ability possessed by 
different kind of light beams, thus extending our previous work on this subject [10]. For 
all the scalar beams studied here, we model the obstacle as a soft-edge aperture with a 
Gaussian profile, in order to be able to produce formulas in closed form. In particular, after 
a careful analysis of the self-healing mechanism from a mathematical and physical point of 
view, we can hnd a new expression for the minimum reconstruction distance Zmin- Then, a 
witness function that quantihes the similarity between the field of the unperturbed beam 
(namely, the beam that would propagate as if the obstacle were not present) and the held of 
the perturbed one (that is, the beam that propagates behind the obstruction), is proposed 
and tested. Finally, our hndings are illustrated by means of a quantitative comparison 
between the self-healing properties of a Gauss beam and a Bessel beam. The notation 
utilized throughout this paper is given in Appendix 


II. SELF-HEALING AND THE EIGENVALUE PROBLEM 
A. Basic definitions and formal developments 

Gonsider a scalar held f{x, y, z) propagating along the 2 ;-axis. An obstruction, character¬ 
ized by a given amplitude transmission function to{x,y), is placed in the plane z = 0. The 
amplitude fo{x, ?/, 0) of the obstructed held in the plane z = 0 can be written as 

fo{x,y,0) =toix,y)f{x,y,0). (1) 

According to the dehnition of angular spectrum, the amplitude fo{x,y, z) of the held trans¬ 
mitted at distance ^ from the obstruction can be written as 

OO 

foix, ^ / / Po{kx, ky) exp [i {xk^ -1- yky -f zk^)] dk^dky, (2) 
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where with = kl + k^ and Fo{k^,ky) = F [f{x,y,0)to{x,y)] {k^,ky), 

namely 


-fo ( kx 1 ky ) 


^11 f{xF,0)to{x,y)exp[-i{xka:+ yky)] dxdy 



— OO 


(3) 


where d'^n = dn^duy, d'^K,' = dF^dn'y, et cetera. The hnal expression for fo{x,y,z) is thns 
given by 


fo{x,y,z) 


(2vr) 


exp {ip - R) exp {izkz 



K')F{K')d^K' 


d^K. 


(4) 


Given the fnnction to{x,y), one can always dehne the transmission fnnction tA{x,y) of an 
apertnre complementary to the obstrnction via the relation (Babinet principle) 

tA{x,y)+to{x,y) = 1. (5) 


By dehnition, both tA{x,y) and to{x,y) are non-negative real-valned fnnctions. Using this 
eqnation into Eq. Q yields to 

fo{x,y,0)= [l-tA{x,y)]f{x,y,0) 

= f{x, y, 0) - tA{x, y)f{x, y, 0) 

= f{x,y,0) - fA{x,y,0). ( 6 ) 


Snbstitnting Eq. ([^ into Eq. Q one obtains 


fo{x,y,z) = 


(27r) 


exp {ip - R) exp {izkz) 


X F{R) - 11 Ta{R - R')F{R') d\' 


d^K 


= f{x,y,z) - fA{x,y,z), 


(7) 
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where Ta{J^) = J^\tA\{kx,ky) and the definition of the “transmitted” fnnction fAix,y,z) is 
obvious. Taking the absolute value squared of both sides of either Eq. ([^ or Eq. Q (the 
integrals are anyway independent from z) and integrating over the whole xj/-plane, we obtain 


I{fo) = I{f - Ia) 

= ^(/)+^(/a) - 2Re II f*{x,y,0)fA{x,y,0)dxdy, 


( 8 ) 


where the average beam intensity I{h) is defined as 


k{h) = / / h*{x,y, z)h{x,y, z) dxdy. 


(9) 


Equation (|^ already shows that perfect self-healing is impossible even in principle because 
the intensity of the transmitted field is unavoidably reduced unless /(/a) = 0. 


B. Defining self-healing 


An optical beam is dubbed “self-reconstructing” when has the ability to recover its initial 
amplitude or intensity profile after interaction with an obstacle. This means that the field 
of a self-reconstructing beam is expected to obey the law 




( 10 ) 


for z > zq, where zq denotes the so-called minimum reconstruction distance and the scaling 
factor Ao = [/(/o)//(/)] accounts for the average intensity reduction caused by the 
interaction with the obstruction. The left side of Eq. (10) is given by Eq. (|^, while the 
right side can be written as 

OO 

Ao/(x, ^ ky) exp [i {xk^ -F yky + zk;,)] dk^^dky 


— OO 
OO 


= ^ / / exp (ip- R) exp {izkf) XoF{k) 


d^K. 


( 11 ) 


Substituting Eqs. (|^ and ( [II| into (10), we obtain the following simple relation: 

Fo{k) k. \qF{k), 

where Eq. ([^ has been used. 


( 12 ) 


It is crucial to notice that Eq. (12) does not contain the variable ^ while, conversely. 


the relation (10) is supposed to be true only for 2 ; > zq. The latter requirement cannot be 


4 








ignored because by definition Eq. (10) cannot be satisfied at z = 0 where, instead, Eq. (§ 
must be fulfilled. Hence, we are apparently faced with a paradox here. In fact, Eq. (12) 
constitutes more a statement about the obstruction rather than the field. This can be seen 
by using Eq. (0 to rewrite ([I^ in the more suggestive form 


1 



— OO 




(13) 


If we replace the approximation symbol with the equality one “=” in Eq. (13), then 
the latter takes the form of a homogeneous Fredholm integral eguation of the second kind 
[m with the unknown function F{k), which has to be an eigenfunction, associated with the 
eigenvalue Aq, of the integral kernel To{k — k') describing the obstruction. This means that 
the requirement (10) is indeed too much restrictive because it can be satisfied only by those 
beams whose angular spectrum (the eigenfunction) is unaffected by the interaction with the 
obstruction, apart from a trivial proportionality factor (the eigenvalue), as shown in Eq. 



FIG. 1. The orthogonal projection of the obstruction on the xy-plane is represented by the 
dark-gray area. This region is circumscribed by the dashed circle of radius a (circumradius) and 
it inscribes the dot-dashed circle of radius b (inradius). Both circles are centered along the z-axis 
of the beam at x = 0 = y. 


C. Redefining the minimum reconstruction distance 

Key to the apparent paradox illustrated above, is the existence of the minimum recon¬ 
struction distance zq after which a self-reconstructing beam restores its initial profile. For 
a single plane wave exp{ik ■ r) with wave vector k = k{x sin 6 cos 4> + y sin 6 sin (p + z cos 6), 
this parameter can be straightforwardly determined in the context of either geometrical 
and wave optics |T0]. Consider an obstruction whose orthogonal projection on the xj/-plane 
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occupies an area O, and let a be the radius (circumradius) of the circle inside which such 
projection can be inscribed and centered on the axis of the beam (see Fig. above). Then 
Zq can be estimated as 


Zo oc 


a 

tan^’ 


(14) 


where the proportionality factor essentially depends on the shape of the obstruction. How¬ 
ever, a beam of light can always be thought as a bunch of plane waves whose density 
in the fc-space is determined by the absolute value squared \F{kx, ky)f of the angular 
spectrum of the held f{x,y,z) representing the beam. Moreover, given the wave vector 
k = k {x sin 9 cos cj) + y sin 9 sin cj) + z cos 9) = xkx + yky + zkz, it is evidently possible to 
rewrite 


1 

tan^ 


kz 

(ki + kiY'" 

(e -ki- kiY'^ 
(ki + kiY'"- 


(15) 


provided that k^ + ky < k‘^. This condition is necessary to maintain kz real-valued and it 
limits the applicability of the equation above to beams whose angular spectrum does not 
contain evanescent waves [12]. From Eq. ( [I5| it follows that we can regard the expression 
of Zo given in Eq. (14) as a function of k = {k^ + kyY^'^ in the fc-space, namely 


Zq —^ a Zi^tY 


K, 


(16) 


On account of the fact that the transverse wave vector k has a density distribution function 
\F{kx, ky)f, we (arbitrarily) define the minimum reconstruction distance Zq as the expected 
value of the function aZ{K), namely 


{kf — K 


2n1/2 


\F{K)f d^K 


K 


Zo irj, 

— = ^ « =- 

a 


|F(k)|^ 


(17) 


ih2 J_ t.2<rJU2 


where both integrals are limited to the disk of equation k‘1 + ky < k'^. It is convenient to 


express the integrals in Eq. (17) in cylindrical coordinates k^ = ncoscp, ky = k sin cp to 
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obtain 


a 


2n 1 

J d^p J ds (^1 — J^{s,ip) 

_0 _ 0 _ 

27r 1 

J dip J ds s J^{s,p) 

0 0 
2tt 7 r /2 

y (i(^ J dO cos^9 {sin 9, p) 

0 _ 0 _ 

2w 1^12 ' 

J dp J d9 sin 9 cos 9 {sin 9, p) 

0 0 


(18) 


where we have defined s = n/k = sin 6* and T{s,p) = |F(kcos 99, Ksin(^)|^. 

In the remainder of this section, we check the validity of Eq. ( [I^ for the cases of i. a 
fundamental Gaussian beam and 2. a Bessel-Gauss beam] both obstrncted by a soft-edge 
Gaussian aperture. For the sake of clarity in the following examples we shall restrict our 
attention to the paraxial regime of propagation. 


1. Fundamental Gaussian beam 


Consider the transmission of a fundamental Gaussian beam of waist Wq across a soft- 
edge Gaussian obstacle of full width 2a located along the axis of the beam ai z = 0. The 
obstruction is described by the transmission function 

to{x,y) = l-exp\- ^^ 1, (19) 


2a2 

where po = xxq + yyo represents the displacement of the obstacle with respect to the beam 
propagation axis. A straightforward calculation gives 

,2 


To{kx, ky) = 27 t6 {k) — exp {—i k ■ po) exp 


-yG+k) 


( 20 ) 


The field describing the Gaussian beam can be written as f{x,y,z) = exp {ikz) g{x,y, z), 
with g{x, y, z) being the fundamental solution of the paraxial wave equation 

1 


9{x,y,z) = 


exp 



fx^ + y^\ 

2 

\z - izRj _ 


( 21 ) 


Z - tZR 

and zr = kw^/2 denotes the Rayleigh range. The Fourier transform at = 0 of this field 
can be easily calculated and the result is 


G{kx, ky) = exp {kl + kl) 


( 22 ) 
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From Eqs. ( 3pl 22) and a straightforward Gaussian integration, we obtain 

k IpoI^ 


^Aikxi ky') 


la 


a^k + zr 


exp 


2 a^k + Zr 


-IK- Po 


Zr 


a^k + Zr 


exp 


a^kl + kl 

y ' (23) 


2 a‘^k + Zr 


with, by dehnition, Go{kx, ky) = G{kx-i ky) — GA{kx, ky). Using this result into Eq. (|§ yields 
to the following expression for the beam transmitted beyond the obstacle: 


9A{x,y,z) = 


Kr 


ZrZ- IttR 


exp 



fx‘^ + y^\ 

[2 

\z - iaRj _ 


(24) 


where, for the sake of clarity, we have chosen po = 0 and we have dehned the modihed 
Rayleigh range kr of the beam gA{x,y,z) transmitted by the aperture complementary to 
the obstruction, as 

< Zr. (25) 

1 + — 
ka^ 

The self-healing capability of a fundamental Gaussian beam is vividly illustrated in Fig. 



FIG. 2. Plots of the intensity distributions (absolute value squared) evaluated at y = 0, of (left to 
right): the incident field g{x,0,z), the “virtual” field transmitted by the aperture complementary 
to the obstruction gA{x, 0, z), and the field transmitted behind the obstacle go{x, 0, z). The plots 
are generated for a Gaussian beam with angular aperture Oq = 2/{kwo) = 7r/12 ~ 15° and a 
soft-edge Gaussian obstruction with full width a/wQ = 0.2 and xq = 0 = yo- At zJzr = 2 the 
intensity profiles of the fields y(x, 0, z) and go{x, 0, z) appear very similar. 


A close inspection of this hgure together with Eqs. (24 25) reveals how the mechanism 
underlying the self-reconstruction process works. From Eq. (25) it follows that qr < zr. 
Ergo, the “virtual” held gA{x,y,z) transmitted by the complementary aperture spreads in 
the xy-plane, while propagating along the z-axis, much more rapidly than the unperturbed 
held g{x,y,z). Therefore, ior z > 2 the intensity prohle of the obstructed beam almost 
coincides with the prohle of the unperturbed one. This process is depicted in Fig. where 























the normalized difference A{z) of the on-axis field intensities /(0,0,z) = | 5 f( 0 , 0, and 
Ja( 0,0,2:) = |5fA(0, 0, is plotted as function of z/zr\ 


M \ -Ia{ 0 , 0 ,z) 

noAz) 


(26) 


From Eq. 


(22) it follows that iF{s,(p) 


|G(kcos(P, fi:sin(p)|^ = exp (—25^/00), where 



FIG. 3. Plots of the the normalized difference A{z) of the on-axis field intensities l{0,0,z) = 
1^(0, 0 , 2 )!^ and /a(0,0,2;) = |(7a(0, 0, z)!^. For long propagation distances z ^ zr the asymptotic 
value is A(oo) = 1 — (aR/zR)'^ = (1 + 4a^)/(l + 2a^)^, where a = ajwQ. 


9o = 2/{kwo) denotes the so-called angular spread of the Gaussian beam [12]. Using this 
result in Eq. (18) yields to 


zq vr i,{i/e^) + h{i/e^) 
a 29 Q sinh(l/ 6 *g^) 

_ (27r)^/^ 1 

6*0 1 - exp (- 2 / 002 )’ 


(^0 < 1), 


(27) 


where lu{z) denotes the modified Bessel function of the first kind of order u [13]. A plot of 
Zo/a is given in Fig. a). Since for 0o ^ 1 if has 0o ~ tan0o, then in the paraxial regime 
of propagation 


zo (27r)^'^^ 

_ _ ___ 

a tan 9q 


(00 ^ 1), 


which is consistent with the expected geometrical optics result. 


( 28 ) 
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b) Bessel-Gauss, in = 0 



FIG. 4. Plots of the minimum reconstruction distance zq/q as given in Eq. (18), evaluated for 


a) a Gaussian beam and h) a zeroth-order (m = 0) Bessel-Gauss beam, a) Plot of Eq. (27) 


(green line) compared with the modified geometrical-optics value (black line), h) Scatter plots of 


numerically integrated Eq. (18) with J-{s,ip) given by Eq. (31) for different angular spreads 9q of 
the Gaussian envelope. Green circles: 9q = 0.5°; blue squares: 9q = 1.0°. The black line represents 
the “one-over-tangent” geometrical-optics prediction and is presented for comparison. 


2. Bessel-Gauss beam 


Consider a Bessel-Gauss beam transmitted across the same soft-edge Gaussian obstruc¬ 
tion as above. The held describing such beam can be written in = 0 as 


/(x,?/, 0) = (±l)™exp (±im0) Jm(Kop)exp ( 


Wr 


(m = 0,l,2,---), (29) 


where p = (x^ -|- kq = /c sin-do and Ju{z) denotes the Bessel function of the hrst kind 

of order u 0. The two key parameters characterizing the beam are the aperture do of the 
Bessel cone in fc-space, and the waist Wq (or, equivalently, the angular spread Oq = 2/{kwo)) 
of the Gaussian envelope. The Fourier transform at z = 0 of this held is given, for m > 0, 
by the following expression: 


w. 


F{k^, ky) = — exp {±imip) I, 


K Ko 


exp 


+ Kq 


4/wl 


Taking the absolute value squared of the function above yields to 


= I; 


2 sin do 

^;r 


exp 


^ (s^ + sin^ i)o) 


e, 


(30) 


(31) 


In this case the integrals in Eq. (18) cannot be calculated analytically and a numerical 


evaluation of the latter is necessary. The resulting values of Zo/a are portrayed in Fig. 1 6) 
as functions of the aperture do, non necessarily paraxial, of the Bessel cone. 
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For 00 ^ 1 the Bessel-Gauss beam reduces to an ordinary Bessel beam. In such a limit 
we recover the expected geometrical optics result: 


a 


tan "00 


(»0 ■« 1 ) ■ 


(32) 


The two examples shown in this section show that Eq. (17) furnishes an appropriate 
measure of the minimum reconstruction distance Zq. 


III. RESOLUTION OF THE PARADOX 


In this section we seek a solution for the paradox outlined earlier. The conundrum may 
be stated as follows: How is it possible to obtain the simultaneous validity of both Eq. ([^ 
and Eq. (10)? In mathematical terms, the problem amounts to understand if and how it is 
possible achieve both 


fo{x,y,0) =to{x,y)f{x,y,0) AND fo{x,y, z) ^ Xof{x,y, z) ^z > Zq. (33) 

As it will be clear soon, the solution of this problem is closely connected to the question 
raised by Chu&Wen [9] about how to quantify the similarity between the two functions 
fo{x,y,z) and f{x,y,z) in order to describe the self-reconstruction ability of a beam. The 
reader is addressed to Appendix for a shortly review of the Chu&Wen approach. 


A. Similarity as relative distance 


Let us begin our discussion about similarity between functions by illustrating a simple, 
prototypical example of a self-reconstructing beam. Consider a “skeleton” Bessel beam made 
of two plane waves only, with wave vectors lying on the xz-plane and forming the angles 
’do and —'do, respectively, with the z-axis. Let f{x,y,z) be the complex-valued scalar field 
representing such a “beam” for z < 0: 


f{x, y,z) = - exp [ikz cos ’do) [exp [ikx sin do) + exp {—ikx sin 'do)] • 


(34) 


This beam hits ai z = 0 a. semitransparent obstacle described by the Gaussian transmission 
function (19) that here we rewrite as 

toix,y) = 1 - exp 


(a;2 - xo) +{y- yoY 


2a2 


(35) 


with a > 0. The held fo{x,y,z) transmitted beyond the obstacle at z > 0 can be straight¬ 
forwardly calculated and the expression is 


fo{x,y,z) 


1 

1 (x - xo)^ + {y- yof 

exp 

. k sin^ 'do z 

1 + iz/ika'^) 

2a? l + iz/{ka?) 

2 l-\-iz/{ka^) 


X cos 


fcsin'do 


X + ixoz/{kad) 

1 -I- izj{ka?) 


(36) 
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It is easy to verify that, as expected, 


lim /o(x, y, 0) = /(x, y, 0). (37) 

0.^00 

The “virtual” held fA{x,y,z) transmitted by the aperture complementary to the obstacle 
can be obtained from the relation 


f{x, y, z) = fo{x, y, z) + Ja^x, y, z). (38) 

The absolute value squared of the three helds /, Ja and fo is shown in Fig. below. 



FIG. 5. Plots of the intensity distributions (absolute value squared) evaluated at y = 0, of (left to 
right): the original field f{x, y, z), the “virtual” field transmitted by the aperture complementary to 
the obstruction fAix, 0, z), and the field transmitted behind the obstruction fo{x, y, z). The plots 
are generated with i9o = 7r/32 ~ 6°, ka = 5 and xq = 0 = yo- For these values the geometrically 
predicted minimum reconstruction distance is equal to zq = tan?9o ~ 72/A:. At kz = 100 the 

intensities of the fields f{x,y,z) and fo{x,y,z) appear already very similar. 


From this hgure it is evident the “self-healing mechanism” in action: during propagation 
from z = 0, where the obstacle is located, to kz = 100 the held practically recovers its orig¬ 
inal intensity distribution. However, it should be noticed that due to the hnite transverse 
extent of the obstruction, the intensity lo{x,0,z) = |/o(x, 0, 2 ;)|^ changes considerably dur¬ 
ing propagation only within the region |A:x| < 30. Conversely, it remains constantly close to 
l{x,0,z) = |/(x,0,2;)|^ in the complementary region |A:x| > 30. This phenomenon appears 
more clearly if one displays in the same hgure the three intensity distributions at diherent 
z, as shown in Fig. [^next page. 

This means that for all practical purposes Eq. (10) represents a far too much restrictive 
constraint. What one really needs is simply to satisfy (10) on the x|/-plane in the neigh¬ 
borhood of the propagation axis z. This statement may be formalized as follows. Consider 
again an obstruction whose orthogonal projection on the x|/-plane occupies the region O, 
and let E be an arbitrary area in the xy-plane strictly contained within O, namely E <Z O. 
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a 


k z = 10 


] kz = 0 b) 




f) kz =20 


d) k z = 40 




FIG. 6. Intensity distributions I{x,y,z) = \f{x,y,z)f (red lines) and Io{x,y,z) = \fo{x,y, z)f 
(green lines) evaluated at y = 0 and a) z = 0, b) kz = 10, c) kz = 20 and d) kz = 40. The black 
dashed curves represent the intensity transmission function \to{x,y)\^ evaluated at y = 0. The 
values of the parameters AQ,i}o,a,xo and yo are the same as in Fig. It is evident how at any 
z > 0 it has Io{x, 0, z) ~ I{x, 0, z) for k |x| > 30 = 6ka. 


For example, E can be the region confined by the inner circle of radius b in Fig. [T} Then, 
as necessary condition for self-haling, we require that the amplitude fo{x,y,z) of the ob¬ 
structed beam is proportional to the amplitude f{x, y, z) of the unperturbed beam only 
within E: 


fo{x,y,z) 


{x,y)&E 


^of{x,y,z] 


{x,y)&E 


'iz > Zq. 


(39) 


From a mathematical point of view, Eq. (39) makes much more sense than Eq. (10). In 


fact, according to the theory of angular spectrum representation of a beam [12], the held 
conhguration at z = 0 completely determines the held distribution at z > 0. Therefore, 


if at a certain distance 2 ; Eq. (10) were satished upon all the xy-plane, then it should be 


also valid at z = 0. But the latter statement is clearly false because at ^ = 0 one has, by 
dehnition. 


/o(i, V, 0) = to{x, v)f(x, y,:) 4 f(^, y, 0- 


(40) 


Thus, we have shown that the origin of the apparent paradox (33) resides in the desideratum 
of satisfying both equations ([^ and (10) over all the xy-plane. This is also the reason why 
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the similarity defined in equation (1) of Ref. [9], fails to furnish a quantitative description 


of self-healing: The double integral defining the scalar product (Bl) extends upon the whole 
a;j/-plane and this erases the 2 ;-dependence. 

To circumvent this difficulty, in this work we propose to define the scalar product in the 
space of functions L 2 {E) as 


(/,^)= / f*{x,y,z)g{x,y,z)dxdy, 

Je 


(41) 


where the integration is now restricted to the domain E. With this definition, the scalar 


product (/, S') naturally becomes a function of 2 ;. Of course the definition (41) is to some 


extent arbitrary in that the choice of the integration domain E is partially discretionary 
(the only constraint is to be entirely contained within O). However, it is useful to remind 
here that the concepts themselves of “self-healing” and “minimum reconstruction distance” 
suffer from the same kind of arbitrariness. In other words, since both Eqs. ([^ and (10) are 
impossible to satisfy over the whole a:|/-plane, one is forced to chose where these equations 
should be satisfied. A reasonable choice is to take E as the region bounded by the inner 
circle of radius 6 < a in Fig. However, different symmetries in the problem may dictate 
different choices, as we will explicitly show later in two examples. 

Motivated by the introduction of the scalar product €D and recalling that the distance 
d{f,g) between two functions in L 2 {E) can be defined as d{f,g) = \\f — g\\, where ||/|| = 
(/./)■'" la , we found it convenient to introduce as z-dependent witness of self-healing, the 
function 


w{z) = 1 — dr, 

where the relative difference dr is defined as 


(42) 


dr = 


11/ -/oil 

+ II/' 


o 


UaZa)''" 


(/./)“'" + [(/./)- 2 Re (/, u) + (/.,, /.4)|'''"' 


(43) 


The limiting values of this witness function w{z) are simply evaluated as follows. For a 
totally opaque obstacle, at 2 ; = 0 it must be /o = 0 for all points {x, y) G E. Therefore, 


from the definition (43) it follows that r(;(0) = 0. Vice versa, if for z > Zq it has fo ~ Aq/ 
for (x, y) G E, then 

/I _ A \ V2 

w{z) 1 — ( --^ 1 for Ao —)■ 1. (44) 


1 +Ac 


Therefore, we have 


0 < w{z) < 1. 


(45) 


14 









In the remainder of this section we will test the effectiveness of our witness function by 
using it to asses the self-healing ability of: 1. the two-plane-wave “beam” described above, 
and 2. a Gaussian beam. 


1. Two-plane-wave beam 


The two-plane-wave held (34) does not depend upon the variable y and has essentially a 


Cartesian geometry. Hence, we choose for E a square of side L centered at (x = 0,?/ = 0), 
namely E := {{x,y) G : |a;| < L A \y\ < L}. In this case the witness function w{z) can 
be calculated analytically. However, the result of this calculation is very cumbersome and 
for the sake of clarity it will not be reported here. In Fig. [^we display w{z) as a function 
of z for different sizes of the square region E. 


k a = 10, 


i9o = ^112 



FIG. 7. Plot of the witness function w{z) given in Eq. (42) evaluated for the two-plane-wave 


field (34) and the Gaussian obstruction (35) of width ka = 10 and centered at xq = 0 = yo- 
The aperture of the “beam” is = 7r/12 and zq ~ 58.5/k denotes the minimum reconstruction 
distance calculates as zq = R{a)/ tan0, where R{a) = a [in (6 -|- 2®/^)] denotes the half-width at 
half maximum of the intensity transmission function \toix,0)\‘^. It is evident that when L increases, 
the dependence of w{z) from z becomes weaker and weaker. The continuous blue line represents 
the asymptotic value of w{z) achieved for L/o —)• 0. 


Since a soft-edge Gaussian obstacle has not a sharp boundary, we can shrink the region E 
to the single point x = 0 = y and obtain the asymptotic value for the witness function 
represented by the continuous blue line in Fig. [^and denoted with w{z). 

The geometrical-optics value for the minimum reconstruction distance given in the liter¬ 
ature, is obtained by calculating the point on the z-axis, namely the point at x = 0 = y, 
where the intensity of the obstructed held begin to raise from zero. Therefore, we can use the 
asymptotic form iv{z), which is indeed calculated “on-axis”, to estimate zq as follows. First, 
we choose a “threshold” value, say w = A, for the witness function. Then, we consider the 
beam as reconstructed only for those distances z from the obstruction such that w{z) > A. 
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By (numerically) inverting this relation we find the minimum reconstruction distance zq as 


zo = w \A). 


The plots of w ^(A) as functions of either i9o 
compared with the geometrical-optics values. 



(46) 

and a are presented in Fig. where they are 


b) do^zrin 



FIG. 8. Plots of '(D“^(A) in Eq. (46) as functions of either a) the aperture and 6) the width of 
the transmission function a. a) Plot of u)“^(A) (blue line) as a function of '&o, with A = 0.8 and 
ka = 20. For comparison, the plot of the geometrical-optics value zq = R{a)/ tan0 (green line), 
is shown, b) Plot of u)“^(A) (blue line) as a function of a, with A = 0.8 and '&o = 'KjVl. For 
comparison, the plot of the geometrical-optics value zq = R{a)l tan0 (green line), is shown. 


In both cases the agreement between w ^(A) and the geometrical-optics value of zq is ex¬ 
cellent. 


2. Gaussian beam 


The Gaussian beam studied in Sec. IICl whose field is given by Eq. (21), possess a 


cylindrical symmetry about the propagation axis z. Therefore, now we choose for E the disk 
of radius h < a depicted in Fig. Also in this case w{z) can be calculated analytically. In 
Fig. 1^) we display w{z) as a function of z for different values of b. When b goes to zero, 
we obtain the simple asymptotic form 


w(C) = 1 - 


(1 + 0 


2 \ 1/2 


C 

2o2 


1 + Of (1 + 2(^2)" 


1/2 ’ 


(47) 
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where ( = zj zr and a = a/ wq. It is interesting to notice that this fnnction is “nniversal” in 
the sense that it does not depend explicitly on the angular spread 6q of the Gaussian beam. 
The plots of w{z) as function oi zj zr for different values of a/tco, are presented in|^). It is 
evident that when the size a of the obstruction goes to zero, the witness function tends to 
the unity. 


a = 0.25 wq b - 0 




FIG. 9. a) Plots of the witness function w{z) given in Eq. (42), evaluated for the Gaussian 


field (21) of waist wg and different radii b of the integration region E. It is evident that when 


b increases, the dependence of w{z) from becomes weaker. The continuous blue line represents 
the asymptotic value •w{z) achieved for 6 —)• 0. b) Plots of the asymptotic witness function w{z) 


(47) for several values of the width a of the soft-edge Gaussian obstruction. When the obstruction 
shrink to zero, then w{z) —)> 1, as expected from physical considerations. 


In order to see the practical signihcance of w{z), in Fig. 10 next page we plot four sections 


of the unperturbed Gaussian beam (21) at different distances z from the obstacle, and we 


compare it with the obstructed beam prohle. The witness function w{z) evidently provide 
for a quantitative estimation of the similarity between these two helds. 


IV. COMPARISON BETWEEN A GAUSSIAN BEAM AND A BESSEL BEAM 


In this last section we compare a Gaussian with a Bessel beam both transmitted behind a 
soft-edge Gaussian obstruction and propagating in the paraxial regime. For comparison, we 
take the cone aperture 'dg of the Bessel beam equal to the angular spread of the Gaussian 
beam. In Figs. 11 and 12 we plot the intensity distributions (absolute value squared) 


evaluated at y = 0, of (left to right): the incident held, the “virtual” held transmitted by the 
aperture complementary to the obstruction, and the held transmitted behind the obstacle. 
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a) 


z - 0, w{z) » 0 


z^zrIIO, w(z)^0.6S7 




x/vi’o x/M'o 


c) z=zrI5, w{z) ^ 0.810 z = zr, w{z) ^ 0.945 




-vAi’o x/m'o 


FIG. 10. Intensity distributions I{x,y,z) = \g{x,y, z)\‘^ (red lines) and Io{x,y,z) = \go{x,y, z)f 
(green lines) evaluated at y = 0 and a) z = 0, b) kz = 10, c) kz = 20 and d) kz = 40. The black 
dashed curves represent the intensity transmission function \to{x,y)\'^ evaluated at y = 0. Here 
^0 = vr/12 and a/wo = 0.2. 


Next, we plot the asymptotic witness functions wg{z) and wb{z) for the Gaussian and the 
Bessel beam, respectively. The witness function is dehned by the values of the intensity on 
the 2 ;-axis according to the formula 


w{z) = 1 — 


V^Aiz) 


+ V^o{z) 


(48) 


where I#{z) = |/^(0,0,^)| . In Fig. 13 one can see that wb{z) > wg{z) always, thus 


indicating a better self-healing property. The vertical line represents the geometrical-optics 
prediction for the minimum self-reconstruction distance zq/zr ~ 0.312, with 


wg{zo) ^ 0.87, wb{,Zo) ^ 0.90. 


(49) 
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|^o(A-/),r)p 



FIG. 11. The plots are generated for a Gaussian beam with angular aperture Bq = 2/{kwo) = 
7r/32 ~ 6° and a soft-edge Gaussian obstruction with full width a/wo = 0.2 and xo = 0 = yo- 



FIG. 12. The plots are generated for a Bessel beam with cone aperture = 7r/32 ~ 6° and a 
soft-edge Gaussian obstruction with full width a/wo = 0.2 and xq = 0 = yo, where wq = 2/{k'&o). 

The two horizontal lines represent the “hdelity” values 0.9 and 0.95 such that 

zoB = Wb^(0.9) ^ 0.298, zoB = iUs\0.95) ^ 0.597, (50) 

zqg = ~ 0-413, zoG = 'u;g^(0-95) ^ 1.20. (51) 

These results show that the Gauss beam reconstruct itself after the Bessel beam. Moreover, 
there is a limit to the reconstruction ability of a Gauss beam, because 

lim wg(,z) = ^ 3 , lim wb{z) = 1, (52) 

Z/ZR^OO 1 -|- y z/zr^OO 

so the hdelity for a Gaussian beam cannot reach the value 1. This is understandable because 
the Gauss beam has a hnite energy which is partially lost because of the obstruction. 
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00 = 7rl32 



ZjZR 

FIG. 13. The plots are generated for Bessel and Gauss beams with 6q = -do = 7r/32 ~ 6° and a 
soft-edge Gaussian obstruction with full width a/wo = 0.2 and xq = 0 = yo, where wq = 2/{k6o). 

V. CONCLUSIONS 

In this work we studied and compared the self-healing properties of Gaussian, Bessel, 
Bessel-Gauss beams and of a pair of plane waves. We found a novel “universal” dehnition 
for the minimum reconstruction distance based on the angular spectrum distribution of a 
beam. Numerous examples are given in the text and illustrated with explanatory hgures. 


Appendix A: Notation 


Three-dimensional vectors in either real and Fourier space are denoted with Latin letters: 
r = XX + yy + zz, k = k^x + kyj) + k^z. Two-dimensional vectors in either real and Fourier 
space are denoted with Greek letters: p = xx + yy, k = k^x + kyj). Gylindrical coordinates in 
real and Fourier spaces are denoted with (p, 0, z) and (k, p, fcz), respectively, with 
and = kl + kf,. All fields considered here are monochromatic with wave number k and 
angular frequency u = ck. 

The Fourier transform of a function g{x, y) of two independent variables x and y will be 
denoted either by J^[g\{kx, ky) or by G{kxi ky) and is dehned by 


G{k^, ky) = — II g{x,y) exp [-i {xk^ -f yky)] dxdy. 


(Al) 


Similarly, the inverse Fourier transform of a function G{kx, ky) will be represented either 


20 








by T ^[G](a;, I/) or by g{x,y) and is defined as 


g{x, y) = ^ II ky) exp [i {xk^ + yky)] dk^dky. 


(A2) 


The two-dimensional Dirac delta symbol d{K — k') stands for 

5{^fl-R')=5{k,-k'Mky-k'y), 

and is defined as 

OO 

1 


(A3) 


5{k — k') = 


(2i) 


exp [ix {kx — k'^)] exp [iy {ky — dxdy. (A4) 


Appendix B: The Chu&Wen proposal 


In the article [9], Chn&Wen defined the scalar prodnct in the space of fnnctions L 2 


as 


if, 9 ) = f* {x, y, z)g{x, y, z) dxdy. (Bl) 


Later, in their eqnation (1), they proposed the following measnre of the similarity between 
two given fnnctions / and g\ 

if, 9 ) 


Similarity ~ cos {f ■ 9 ) = 


(B2) 


where the norm of a fnnction is defined as ||/|| = This notation may seems not 

very appropriate becanse the right side of Eq. ( |B2[ ) can be a complex nnmber. A more 
snitable definition yielding to a non-negative nnmber is simply 


Similarity 


(B3) 


which coincides with the standard definition of fidelity of quantum pure states in qnantnm 
mechanics na. 

Using the Parseval’s theorem, Chn&Wen were able to show that the similarity defined as 


in Eq. (B2) cannot depend npon the propagation distance (in their derivation Chn&Wen 


implicitly assnmed that the considered angnlar spectrnm did not contain evanescent waves). 
Therefore, they adopted a new description of similarity, given in their eqnation (11), defined 
as 

(l/l,l^7l) 


Similarity ~ cos (I/I ■ 1^1) = 


ll/ll 11^711 


(B4) 
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With this dehnition, the similarity between / and g becomes a function of the propagation 
distance z. However, if one applies Eq. (B4) to pure (namely, non Gaussian) Bessel beams, 
this similarity function becomes 2 ;-in dependent. Therefore, also the dehnition (B4) in some 
circumstances may be not fully satisfactory. 
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